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ABSTRACT 

The growth of supermassive black holes by merging and accretion in hierarchical mod- 
els of galaxy formation is studied by means of Monte Carlo simulations. A tight linear 
relation between masses of black holes and masses of bulges arises if the mass accreted 
by supermassive black holes scales linearly with the mass forming stars and if the 
redshift evolution of mass accretion tracks closely that of star formation. Differences 
in redshift evolution between black hole accretion and star formation introduce a con- 
siderable scatter in this relation. A non-linear relation between black hole accretion 
and star formation results in a non-linear relation between masses of remnant black 
holes and masses of bulges. The relation of black hole mass to bulge luminosity ob- 
served in nearby galaxies and its scatter are reproduced reasonably well by models in 
which black hole accretion and star formation are linearly related but do not track 
each other in redshift. This suggests that a common mechanism determines the ef- 
ficiency for black hole accretion and the efficiency for star formation, especially for 
bright bulges. 
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1 INTRODUCTION 



, Lynden-Bell (1969) suggested that active galaxies are fu- 
' elled by accretion of matter onto supermassive black holes 
and that massive dark compact objects could therefore be 
, present in the nuclei of many quiescent galaxies as remnants 
' of QSOs. Now there is solid observational evidence for the 
presence of supermassive black holes in our own and in many 
nearby galaxies. A correlation between the mass of a super- 
massive black hole and the luminosity of the bulge of the 
host galaxy was already apparent in the first small sample 
of detections. The larger samples meanwhile available have 
backed up the existence of such a correlation (Kormendy 
1993; Kormendy & Richstone 1995; Magorrian et al. 1998; 
Ford et al. 1998; van der Marel 1998; Ho 1998; Mc Leod 
1998). 

The growth history of these supermassive black holes 
is, however, rather uncertain. It is widely believed that they 
built up during the QSO epoch, around z ~ 2 — 3. As 
pointed out by Phinney (1997) and discussed in more de- 
tail by Haehneh, Natarajan & Rees (1998, HNR98), the 
total integrated mass density in remnant black holes sub- 
stantially exceeds the integrated mass density accreted onto 
supermassive black holes by optically bright QSOs unless 
either the efficiency for producing radiation is rather small 
or the masses of black holes in nearby galaxies are overes- 



timated (see also Richstone et al. 1998 for a review). This 
leaves some room for accretion other than that inferred from 
optically bright QSOs (see also Fabian & Iwasawa 1999 and 
Salucciet al. 1999). 

A striking feature of the activity of optically bright 
QSOs is the fast evolution of their space density. A rapid rise 
between 2 ~ 5 and z ~ 3 is followed by a strong decline after 
2; '^ 2 (e.g. Shaver et al. 1996). This rapid rise can be plau- 
sibly explained by the evolution of the dark matter (DM) 
distribution in hierarchical cosmogonies (Haehnelt & Rees 
1993, HR93): deeper potential wells form at later times and 
the first potential wells which are deep enough to retain their 
gas in spite of the input of energy and momentum due to star 
formation offer the best conditions for efficient formation of 
supermassive black holes. As discussed by HR93 and Cava- 
liere, Perri & Vittorini (1997), the decline in the activity of 
galactic nuclei is probably due to a transition from a high- 
redshift phase, in which supermassive black holes form in ini- 
tially gas rich nuclei, to a low-redshift phase, where QSOs 
are refuelled by mergers of increasingly gas-poor galaxies. 
The late merging of gas-poor galaxies containing supermas- 
sive black holes will obviously affect the mass distribution of 
remnant black holes. For instance, it will introduce scatter 
in any relation between properties of galaxies and masses 
of remnant black holes. Here we use Monte Carlo realiza- 
tions of the merging histories of DM haloes, together with 
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a simple scheme for galaxy formation, to model the merg- 
ing of galaxies in a hierarchical cosmogony. We investigate 
different schemes for the accretion of matter onto super- 
massive black holes and we study the effect on the mass 
distribution of remnant black holes in low-redshift galaxies. 
The structure of the paper is as follows. In Section 2 we 
describe the Monte Carlo simulations. In Section 3 we com- 
pare the resulting relations between masses of black holes 
and luminosities of bulges to that observed in nearby galax- 
ies. Section 4 contains our conclusions. We assume a flat 
cosmology with Qm ~ 0.3, Oa — 0.7 and a Hubble constant 
of 75kms"^Mpc"\ 



2 MONTE CARLO SIMULATIONS OF THE 
MERGING AND ACCRETION HISTORY OF 
SUPERMASSIVE BLACK HOLES 

2.1 Merging histories of dark matter haloes 

The standard paradigm for the formation of structure in the 
Universe is the gravitational instability of density fluctua- 
tions in a primordial Gaussian random field. In this picture, 
the number density of collapsed DM haloes is accurately 
described by the Press-Schechter formula 



N{M,z)dM = 
(l)(27)^(^)(7^)exp(-- 



(1) 



dM 



(Press & Schechter 1974). Here p(z) is the cosmological den- 
sity at redshift z, cro(M)^ is the variance of the linearly 
extrapolated density field on the scale M, a; is defined as 
u}{z) = ScD{0)/D{z), where Sc is the over-density threshold 
at which density fluctuations collapse, and D{z) is the lin- 
ear growth factor of density fluctuations. The variance of the 
density field is related to the power spectrum of the density 
fluctuations. Here we assume a standard cold dark matter 
power spectrum as given by Bond & Efstathiou (1984), with 
a normalisation that reproduces the present-day space den- 
sity of clusters of galaxies (Eke, Cole & Frenk 1996). The 
dependence of the linear growth factor on redshift is given 
by Heath (1977). 

Techniques for generating Monte Carlo realizations of 
the merging histories of DM haloes based on extensions of 
the Press-Schechter formula have been described for instance 
by Cole & Kaiser (1988), Cole (1991), Lacey & Cole (1993), 
Kauffmann & White (1993) and Somerville & Kolatt (1997). 
The basic equation is the conditional probability that a halo 
of mass M at redshift z has a progenitor of mass M — AM 
at redshift z + Az, 



P{M -^ M - AM, z^ z + Az)^ 



(2) 



1 

/2n 



ui{z + Az)—Lj(z) 



[<T2(M-AM)-a2(M)l: 



-3- exp 



l{c^(z + Az)-i^(z)]^ 
'2[<T^(M-AM)-<T^(M)] 



(Lacey & Cole 1993). We follow the description of Somerville 
& Kolatt (1997) to construct Monte Carlo realizations of the 
merging histories of DM haloes from equation (2). The prob- 
ability distribution (2) is used to assign a mass M to a pro- 
genitor at redshift Az of a DM halo of mass Mo at 2: = 0. The 
procedure is iterated and further progenitors are drawn from 



the distribution (2), but the merging history is followed only 
for haloes with circular velocities Vc — (GM/rvir)^ above 
a certain threshold vi. Haloes with Vc < vi are treated as 
accreted mass. Progenitors with a mass larger than the not 
yet allocated mass are discarded. We construct progenitors 
of progenitors until all haloes have Vc < «i . The procedure is 
repeated for a representative sample of haloes at z = 0. We 
adopt a step of Az — 0.01 and a resolution of wi = 70kms^^. 



2.2 The galaxy formation scheme 

The modelling of galaxies within hierarchal cosmogonies has 
reached a considerable level of complexity. Here we concen- 
trate on the effects of mergers on the growth of supermassive 
black holes and on the distribution of remnant black holes. 
Therefore we adopt a rather simple scheme for galaxy for- 
mation similar to that proposed by White and Rees (1978), 
which should nevertheless catch the essential features of the 
hierarchical merging of galaxies (White 1996). 

Important conditions for the formation of a galaxy in 
a collapsed DM halo are the ability of the gas to cool and 
the ability of the halo to retain its gas in spite of the in- 
put of energy and momentum due to star formation and 
supernova explosions. The ability of the gas to cool depends 
on the temperature and the density of the gas. For cooling 
by Bremsstrahlung these dependencies conspire to give an 
upper limit for the mass of the gas which is able to cool 
efficiently (Silk 1977; Rees & Ostriker 1977). Feedback from 
star formation is more important in haloes with shallower 
potential wells and therefore smaller circular velocities (see 
e.g. Kaufltmann, White & Guiderdoni 1993, KWG93). Here 
we do not treat the cooling and feedback explicitly. Instead, 
we introduce an effective efficiency for star formation, which 
depends on the halo circular velocity in such a way that the 
mass of the galaxy associated with a DM halo of mass Mhaio 
at redshift z scales as 



M*(Mhalo) = e*MhaloWc exp [-(Wc/Umax) ] 



(3) 



independently of redshift; e* and Wmax are chosen in such a 
way that the luminosity function of our simulated galaxies 
reproduces the observed luminosity function of galaxies at 
redshift zero (Fig.l). The cut-off in virial velocities mimics 
the inability of the gas to cool and form stars in the very 
deep potential wells which form at late times. 

When DM haloes merge, it will take some time before 
the corresponding galaxies sink to the centre of the merged 
halo due to dynamical friction. To model this, we follow 
KWG93 and assume that initially a single galaxy forms at 
the centre of each halo. When haloes merge, the central 
galaxy of the largest progenitor halo becomes the central 
galaxy of the new halo. The "satellite" galaxies are assumed 
to merge with the central galaxy on the dynamical friction 
time- scale 



^df = 



1.17<i,Wc 



ln(M/Msat)GMs, 



(4) 



(Binney & Tremaine 1987). Here M, rvir, Vc are the mass, 
the virial radius and the circular velocity of the new halo 
and Msat is the total mass of the satellite including its DM 
halo. The stellar mass of the merged galaxy is assumed to be 
the maximum of the sum of the stellar masses of the merging 
galaxies and the mass given by equation (3) . If the latter is 
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Figure 1. Observed and simulated B-band luminosity func- 
tion of galaxies and bulges (data from Salucci et al. 1998 and 
Marzke et al. 1998). A relation of stellar mass to circular ve- 
locity and mass of the DM halo of the form Mi,(M,z) = 
0.014Mhaio X (1)0/200 km s"l)2exp[-(t>c/300 km s"!)'*] is as- 
sumed. Major mergers result in the formation of a bulge for mass 
ratios above /o = 0.2. 



larger, we interpret the difference AM* as tfie amount of gas 
that has formed stars in the disks of the galaxies since the 
last merger. 

We follow KWG93 again and introduce a critical mass 
ratio /e above which a merger produces an elliptical galaxy. 
For mass ratios below the threshold the smaller galaxy 
is assumed to be tidally disrupted and its mass is added 
to the disk of the larger one. The value of /c is deter- 
mined by requiring the Monte Carlo simulations to repro- 
duce the observed bulge luminosity function, once the other 
free parameters have been fixed. A mass-to-light ratio of 
M^,/Lb = 5(Ls/5.8 X IO^'Lq)"-^^ is used (Salucci et al. 
1999). Figure 1 compares the resulting luminosity functions 
to those observed. 



2.3 The growth of supermassive black holes by 
accretion and merging 

In a hierarchical cosmogony supermassive black holes grow 
by merging and by accretion of gas. Here we assume that 
a major galactic merger leads to an immediate coalescence 
of the supermassive black holes contained in the merging 
galaxies. This should be a reasonable assumption for violent 
mergers, though in a smooth stellar background it may take 
longer for the holes to coalesce (Begelman, Blandford & Rees 
1980). Some fraction of the gas which is able to cool in DM 
haloes is assumed to be accreted onto the merged supermas- 
sive black hole. The fraction of cool gas ending up in the hole 
should have a complicated dependence on halo and galactic 
properties (redshift, past accretion history, past star forma- 
tion history, past merging history, angular momentum). We 



Figure 2. The redshift dependence of the comoving mass den- 
sity in supermassive black holes for: (i) Mace = 6 X 10~^AM*, 
(dashed line); (ii) Mace = 1.4x 10~^(1-|-2)2AM* (solid line); (iii) 
Mace = 10"6(l-|-z)2Mhaioexp[-(i;c/300kms"l)'l], (dotted line). 
The points shovif the observed redshift dependence as determined 
by Choski & Turner (1992). 



shall not attempt to model these processes here in detail. 
We shall rather investigate three very simple models for the 
efficiency of accretion. 

(i) During a major merger a fixed fraction of the gas that 
has formed stars since the last merger in the disks of the 
colliding galaxies is accreted onto the supermassive black 
hole formed by the coalescence of the supermassive black 
holes in the merging galaxies (Mace = taecAM*, eace = 6 x 
10-^). 

(ii) As in model (i) but with an accretion efficiency de- 
pendent on redshift to approximate the redshift evolution 
of the accretion onto black holes inferred from the counts 

of optically bright QSOs (Mace = eaec(2)AM*, eace(z) = 

1.4 X 10^^(1 + 2)^). 

(iii) During a major merger a fraction of the total baryonic 
mass of the merged DM halo is accreted onto the super- 
massive black hole with an accretion efficiency dependent 
on redshift as in model (ii) and with the same cut-off for 
high circular velocity as in equation (3) to mimic the in- 
ability of the gas to cool in deep potential wells which form 
at late times (Mace = eacc(2)Mhaio exp [— («c/300kms"^)''], 

eace(2) = 10-«(l+2)2). 

Figure 2 shows the fraction of the total mass accreted 
before redshift z in the Monte Carlo simulations and the ob- 
served fraction inferred from the blue light emitted by op- 
tically bright QSOs (Sofian 1982; Choski & Turner 1992). 
Models (ii) and (iii) approximate the observed redshift evo- 
lution by construction, whereas in model (i) most of the 
mass is accreted at lower redshifts than those indicated by 
the epoch of optically bright QSOs. Model (iii) is similar to 
that recently proposed by Haiman & Loeb (1998). 
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Figure 3. (a) Observed black hole masses and bulge luminosities for the samples of nearby galaxies compiled by Ho (1998, circles) and 
Magorrian et al. (1998, squares). The solid line shows the linear least square fit to the data, log{M g u / Mq) = 1.281og(La/L0) — 4.46, 
while the dashed lines show the itlcT deviation (<t = 0.74). The dashed lines are the same in all panels for comparison, (b) Monte 
Carlo simulations of black hole masses and bulge luminosities with Mace = 6 X 10~^AM<- (model i). (c) Same as (b) with Mace = 
1.4 X 10"-^(1 + z)2AM* (model ii). (d) Same as (b) with Mace = 10"''(1 + z)2Mhaioexp [-(vc/SOOkms^^)-*] (model iii). The linear least 
square fit to the results of model (iii) has a slope of 0.6, shown by the two dotted lines. 



The normalisation eacc is determined by requiring the 
Monte Carlo simulations to reproduce the comoving mass 
density of supermassive black holes at the present time. 
As discussed in the introduction and by HNR98, this 
value is still uncertain, and here we have used a value of 
IO^MqMpc^^. This value is about six times larger than that 
estimated by Choski & Turner (1992) using optical counts 
of unobscured QSOs. 



3 THE MASS DISTRIBUTION OF 

SUPERMASSIVE BLACK HOLES IN 
LOW-REDSHIFT GALAXIES 

Figure 3a shows the observed black hole masses and bulge 
luminosities for a sample of 39 galaxies as compiled by Ho 
(1998) and Magorrian (1998), together with a linear least 
square fit and the ±la deviation. Figures 3b, 3c and 3d 
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show the results of the Monte Carlo simulations. Models 
(i) and (ii) both reproduce the observed roughly linear re- 
lation between black hole mass and bulge luminosity. This 
is mainly due to our assumption that the accreted mass is 
proportional to the mass that has formed stars, which is 
strictly true in model (i) and is still a good approximation 
for model (ii). In model (i) the relation is very tight, tighter 
than observed. In model (ii), instead, the different redshift 
evolution of the mass forming stars and of the mass accreted 
onto supermassive black holes acts as a substantial source 
of scatter, which is now comparable to that in the observed 
relation. Model (iii) produces a shallower relation of black 
hole mass to bulge luminosity; the predicted remnant black 
hole masses in large bulges are here considerably smaller 
than observed; this is due to the non-linear relation between 
bulge mass and halo mass: as the assumed star formation ef- 
ficiency depends on the depth of the potential well, the bulge 
mass rises faster than the halo mass with increasing bulge 
luminosity. The scatter in model (iii) is again comparable to 
the observed scatter and to that in model (ii). 

We have checked the sensitivity of our results to some 
of our assumptions. The slope of the relation of black hole 
mass to bulge luminosity depends only weakly on the cos- 
mology, on the details of the galaxy formation scheme and 
on the quasar epoch once the parameters have been fixed to 
reproduce the present-day bulge luminosity function and the 
integrated mass density in black holes. The scatter in this 
relation, however, grows when the difference between the 
redshift evolution of star formation and black hole accretion 
increases. 



4 DISCUSSION AND CONCLUSIONS 

We have used Monte Carlo realizations of the extended 
Press-Schechter formalism together with simple schemes for 
galaxy formation and mass accretion onto supermassive 
black holes to investigate how the merging and accretion 
history of supermassive black holes within a hierarchical cos- 
mogony affects the mass distribution of remnant black holes 
in low-redshift galaxies. 

The observed roughly linear relation between black hole 
masses and bulge masses can be reproduced if the mass ac- 
creted during a major merger is a fixed fraction of the mass 
of the gas that has formed stars in the merging galaxies since 
the last major merger. This relation is preserved if this frac- 
tion is independent of mass but varies with redshift, but 
such a variation introduces significant scatter into the rela- 
tion. The redshift variation necessary to match the evolution 
of the global mass accretion rate inferred from the integrated 
light emission of optically bright QSOs (assuming that their 
radiation efficiency does not depend on redshift) results in a 
scatter that is comparable to that in the observed relation. 

The linear relation between bulge and black hole mass 
and the modest scatter in this relation at the bright end of 
the bulge luminosity function seem to suggest that for most 
of the integrated mass density in supermassive black holes 
a common mechanism determines the efficiency for mass ac- 
cretion of supermassive black holes and the efficiency for 
star formation. 

The scatter in the observed relation must be at least 
partially due to observational errors and there seems not 



much room for additional sources of intrinsic scatter at the 
bright end of the luminosity function (except in the some- 
what artificial model (i)). At the faint end the scatter in 
the observed relation seems to increase. This might indicate 
that here other parameters, e.g. the spin of the DM halo, 
influence black hole and star formation in a different way. 

If the accreted mass scaled with the mass of the halo 
rather than with the stellar mass, the black hole mass would 
depend on the bulge mass less steeply than observed and 
black holes in luminous bulges would be less massive than 
the observations suggest. 
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